Chaotic escape of impurities and sticky orbits in toroidal plasmas.

We investigate chaotic impurity transport in toroidal fusion plasmas (tokamaks) from the point of view of passive advection of charged particles due to E×B drift motion. We use realistic tokamak profiles for electric and magnetic fields as well as toroidal rotation effects, and consider also the effects of electrostatic fluctuations due to drift instabilities on particle motion. A time-dependent one degree-of-freedom Hamiltonian system is obtained and numerically investigated through a symplectic map in a Poincaré surface of section. We show that the chaotic transport in the outer plasma region is influenced by fractal structures that are described in topological and metric point of views. Moreover, the existence of a hierarchical structure of islands-around-islands, where the particles experience the stickiness effect, is demonstrated using a recurrence-based approach.

Isochronous island bifurcations driven by resonant magnetic perturbations in tokamaks.

Recent evidence shows that heteroclinic bifurcations in magnetic islands may be caused by the amplitude variation of resonant magnetic perturbations in tokamaks. To investigate the onset of these bifurcations, we consider a large aspect ratio tokamak with an ergodic limiter composed of two pairs of rings that create external primary perturbations with two sets of wave numbers. An individual pair produces hyperbolic and elliptic periodic points, and its associated islands, that are consistent with the Poincaré-Birkhoff fixed-point theorem. However, for two pairs producing external perturbations resonant on the same rational surface, we show that different configurations of isochronous island chains may appear on phase space according to the amplitude of the electric currents in each pair of the ergodic limiter. When one of the electric currents increases, isochronous bifurcations take place and new islands are created with the same winding number as the preceding islands. We present examples of bifurcation sequences displaying (a) direct transitions from the island chain configuration generated by one of the pairs to the configuration produced by the other pair, and (b) transitions with intermediate configurations produced by the limiter pairs coupling. Furthermore, we identify shearless bifurcations inside some isochronous islands, originating nonmonotonic local winding number profiles with associated shearless invariant curves.

Fractal and Wada escape basins in the chaotic particle drift motion in tokamaks with electrostatic fluctuations.

The E×B drift motion of particles in tokamaks provides valuable information on the turbulence-driven anomalous transport. One of the characteristic features of the drift motion dynamics is the presence of chaotic orbits for which the guiding center can experience large-scale drifts. If one or more exits are placed so that they intercept chaotic orbits, the corresponding escape basins structure is complicated and, indeed, exhibits fractal structures. We investigate those structures through a number of numerical diagnostics, tailored to quantify the final-state uncertainty related to the fractal escape basins. We estimate the escape basin boundary dimension through the uncertainty exponent method and quantify final-state uncertainty by the basin entropy and the basin boundary entropy. Finally, we recall the Wada property for the case of three or more escape basins. This property is verified both qualitatively and quantitatively using a grid approach.

Fundamental solution of diffusion equation for Kappa gas: Diffusion length for suprathermal electrons in solar wind.

A recent numerical treatment of data obtained by the Parker Solar Probe spacecraft describes the electron concentration in solar wind as a function of the heliocentric distance based on a Kappa distribution with spectral index κ=5. In this work, we derive and, subsequently, solve an entirely different class of nonlinear partial differential equations describing the one-dimensional diffusion of a suprathermal gas. The theory is applied to describe the aforementioned data and we find a spectral index κ≳1.5 providing the widely acknowledged identification of Kappa electrons in solar wind. We also find that suprathermal effects increase the length scale of classical diffusion by one order of magnitude. Such a result does not depend on the microscopic details of the diffusion coefficient since our theory is based on a macroscopic formulation. Forthcoming extensions of our theory by including magnetic fields and relating our formulation to nonextensive statistics are briefly addressed.

Chaotic saddles and interior crises in a dissipative nontwist system.

We consider a dissipative version of the standard nontwist map. Nontwist systems present a robust transport barrier, called the shearless curve, that becomes the shearless attractor when dissipation is introduced. This attractor can be regular or chaotic depending on the control parameters. Chaotic attractors can undergo sudden and qualitative changes as a parameter is varied. These changes are called crises, and at an interior crisis the attractor suddenly expands. Chaotic saddles are nonattracting chaotic sets that play a fundamental role in the dynamics of nonlinear systems; they are responsible for chaotic transients, fractal basin boundaries, and chaotic scattering, and they mediate interior crises. In this work we discuss the creation of chaotic saddles in a dissipative nontwist system and the interior crises they generate. We show how the presence of two saddles increases the transient times and we analyze the phenomenon of crisis induced intermittency.

Large coefficient of variation of inter-spike intervals induced by noise current in the resonate-and-fire model neuron.

Neuronal spike variability is a statistical property associated with the noise environment. Considering a linearised Hodgkin-Huxley model, we investigate how large spike variability can be induced in a typical stellate cell when submitted to constant and noise current amplitudes. For low noise current, we observe only periodic firing (active) or silence activities. For intermediate noise values, in addition to only active or inactive periods, we also identify a single transition from an initial spike-train (active) to silence dynamics over time, where the spike variability is low. However, for high noise current, we find intermittent active and silence periods with different values. The spike intervals during active and silent states follow the exponential distribution, which is similar to the Poisson process. For non-maximal noise current, we observe the highest values of inter-spike variability. Our results suggest sub-threshold oscillations as a possible mechanism for the appearance of high spike variability in a single neuron due to noise currents.

Low-dimensional chaos in the single wave model for self-consistent wave-particle Hamiltonian.

We analyze nonlinear aspects of the self-consistent wave-particle interaction using Hamiltonian dynamics in the single wave model, where the wave is modified due to the particle dynamics. This interaction plays an important role in the emergence of plasma instabilities and turbulence. The simplest case, where one particle (N=1) is coupled with one wave (M=1), is completely integrable, and the nonlinear effects reduce to the wave potential pulsating while the particle either remains trapped or circulates forever. On increasing the number of particles ( N=2, M=1), integrability is lost and chaos develops. Our analyses identify the two standard ways for chaos to appear and grow (the homoclinic tangle born from a separatrix, and the resonance overlap near an elliptic fixed point). Moreover, a strong form of chaos occurs when the energy is high enough for the wave amplitude to vanish occasionally.

Transport of blood particles: Chaotic advection even in a healthy scenario.

We study the advection of blood particles in the carotid bifurcation, a site that is prone to plaque development. Previously, it has been shown that chaotic advection can take place in blood flows with diseases. Here, we show that even in a healthy scenario, chaotic advection can take place. To understand how the particle dynamics is affected by the emergence and growth of a plaque, we study the carotid bifurcation in three cases: a healthy bifurcation, a bifurcation with a mild stenosis, and the another with a severe stenosis. The result is non-intuitive: there is less chaos for the mild stenosis case even when compared to the healthy, non-stenosed, bifurcation. This happens because the partial obstruction of the mild stenosis generates a symmetry in the flow that does not exist for the healthy condition. For the severe stenosis, there is more irregular motion and more particle trapping as expected.

Fractal structures in the parameter space of nontwist area-preserving maps.

Fractal structures are very common in the phase space of nonlinear dynamical systems, both dissipative and conservative, and can be related to the final state uncertainty with respect to small perturbations on initial conditions. Fractal structures may also appear in the parameter space, since parameter values are always known up to some uncertainty. This problem, however, has received less attention, and only for dissipative systems. In this work we investigate fractal structures in the parameter space of two conservative dynamical systems: the standard nontwist map and the quartic nontwist map. For both maps there is a shearless invariant curve in the phase space that acts as a transport barrier separating chaotic orbits. Depending on the values of the system parameters this barrier can break up. In the corresponding parameter space the set of parameter values leading to barrier breakup is separated from the set not leading to breakup by a curve whose properties are investigated in this work, using tools as the uncertainty exponent and basin entropies. We conclude that this frontier in parameter space is a complicated curve exhibiting both smooth and fractal properties, that are characterized using the uncertainty dimension and basin and basin boundary entropies.

Recurrence quantification analysis for the identification of burst phase synchronisation.

In this work, we apply the spatial recurrence quantification analysis (RQA) to identify chaotic burst phase synchronisation in networks. We consider one neural network with small-world topology and another one composed of small-world subnetworks. The neuron dynamics is described by the Rulkov map, which is a two-dimensional map that has been used to model chaotic bursting neurons. We show that with the use of spatial RQA, it is possible to identify groups of synchronised neurons and determine their size. For the single network, we obtain an analytical expression for the spatial recurrence rate using a Gaussian approximation. In clustered networks, the spatial RQA allows the identification of phase synchronisation among neurons within and between the subnetworks. Our results imply that RQA can serve as a useful tool for studying phase synchronisation even in networks of networks.

Riddling: Chimera's dilemma.

We investigate the basin of attraction properties and its boundaries for chimera states in a circulant network of Hénon maps. It is known that coexisting basins of attraction lead to a hysteretic behaviour in the diagrams of the density of states as a function of a varying parameter. Chimera states, for which coherent and incoherent domains occur simultaneously, emerge as a consequence of the coexistence of basin of attractions for each state. Consequently, the distribution of chimera states can remain invariant by a parameter change, and it can also suffer subtle changes when one of the basins ceases to exist. A similar phenomenon is observed when perturbations are applied in the initial conditions. By means of the uncertainty exponent, we characterise the basin boundaries between the coherent and chimera states, and between the incoherent and chimera states. This way, we show that the density of chimera states can be not only moderately sensitive but also highly sensitive to initial conditions. This chimera's dilemma is a consequence of the fractal and riddled nature of the basin boundaries.

Synchronous behaviour in network model based on human cortico-cortical connections.

OBJECTIVE: We consider a network topology according to the cortico-cortical connection network of the human brain, where each cortical area is composed of a random network of adaptive exponential integrate-and-fire neurons. APPROACH: Depending on the parameters, this neuron model can exhibit spike or burst patterns. As a diagnostic tool to identify spike and burst patterns we utilise the coefficient of variation of the neuronal inter-spike interval. MAIN RESULTS: In our neuronal network, we verify the existence of spike and burst synchronisation in different cortical areas. SIGNIFICANCE: Our simulations show that the network arrangement, i.e. its rich-club organisation, plays an important role in the transition of the areas from desynchronous to synchronous behaviours.

Inference of topology and the nature of synapses, and the flow of information in neuronal networks.

The characterization of neuronal connectivity is one of the most important matters in neuroscience. In this work, we show that a recently proposed informational quantity, the causal mutual information, employed with an appropriate methodology, can be used not only to correctly infer the direction of the underlying physical synapses, but also to identify their excitatory or inhibitory nature, considering easy to handle and measure bivariate time series. The success of our approach relies on a surprising property found in neuronal networks by which nonadjacent neurons do "understand" each other (positive mutual information), however, this exchange of information is not capable of causing effect (zero transfer entropy). Remarkably, inhibitory connections, responsible for enhancing synchronization, transfer more information than excitatory connections, known to enhance entropy in the network. We also demonstrate that our methodology can be used to correctly infer directionality of synapses even in the presence of dynamic and observational Gaussian noise, and is also successful in providing the effective directionality of intermodular connectivity, when only mean fields can be measured.

Synchronised firing patterns in a random network of adaptive exponential integrate-and-fire neuron model.

We have studied neuronal synchronisation in a random network of adaptive exponential integrate-and-fire neurons. We study how spiking or bursting synchronous behaviour appears as a function of the coupling strength and the probability of connections, by constructing parameter spaces that identify these synchronous behaviours from measurements of the inter-spike interval and the calculation of the order parameter. Moreover, we verify the robustness of synchronisation by applying an external perturbation to each neuron. The simulations show that bursting synchronisation is more robust than spike synchronisation.

Spike timing-dependent plasticity induces non-trivial topology in the brain.

We study the capacity of Hodgkin-Huxley neuron in a network to change temporarily or permanently their connections and behavior, the so called spike timing-dependent plasticity (STDP), as a function of their synchronous behavior. We consider STDP of excitatory and inhibitory synapses driven by Hebbian rules. We show that the final state of networks evolved by a STDP depend on the initial network configuration. Specifically, an initial all-to-all topology evolves to a complex topology. Moreover, external perturbations can induce co-existence of clusters, those whose neurons are synchronous and those whose neurons are desynchronous. This work reveals that STDP based on Hebbian rules leads to a change in the direction of the synapses between high and low frequency neurons, and therefore, Hebbian learning can be explained in terms of preferential attachment between these two diverse communities of neurons, those with low-frequency spiking neurons, and those with higher-frequency spiking neurons.

Chaotic particle heating due to an obliquely propagating wave in a magnetized plasma.

We study the dynamics of a relativistic charged particle in the presence of a uniform magnetic field and a stationary electrostatic wave that propagates at an arbitrary angle. The wave is considered as a series of periodic pulses which allows us to derive an exact map for the system. In particular, we investigate the heating process of an initially low-energy particle. It is found that abrupt changes in the maximum energy attained by the particle may occur as the angle between the wave propagation and the magnetic field varies. To determine what is the mechanism behind this phenomenon a reduced Hamiltonian that retains the important dynamical features is obtained. Using both Poincaré plots and perturbation theory, we identify that a separatrix reconnection is the key mechanism for the abrupt change in particle response.

Alternate islands of multiple isochronous chains in wave-particle interactions.

We analyze the dynamics of a relativistic particle moving in a uniform magnetic field and perturbed by a standing electrostatic wave. We show that a pulsed wave produces an infinite number of perturbative terms with the same winding number, which may generate islands in the same region of phase space. As a consequence, the number of isochronous island chains varies as a function of the wave parameters. We observe that in all the resonances, the number of chains is related to the amplitude of the various resonant terms. We determine analytically the position of the periodic points and the number of island chains as a function of the wave number and wave period. Such information is very important when one is concerned with regular particle acceleration, since it is necessary to adjust the initial conditions of the particle to obtain the maximum acceleration.

Secondary nontwist phenomena in area-preserving maps.

Phenomena as reconnection scenarios, periodic-orbit collisions, and primary shearless tori have been recognized as features of nontwist maps. Recently, these phenomena and secondary shearless tori were analytically predicted for generic maps in the neighborhood of the tripling bifurcation of an elliptic fixed point. In this paper, we apply a numerical procedure to find internal rotation number profiles that highlight the creation of periodic orbits within islands of stability by a saddle-center bifurcation that emerges out a secondary shearless torus. In addition to the analytical predictions, our numerical procedure applied to the twist and nontwist standard maps reveals that the atypical secondary shearless torus occurs not only near a tripling bifurcation of the fixed point but also near a quadrupling bifurcation.

Effective transport barriers in nontwist systems.

In fluids and plasmas with zonal flow reversed shear, a peculiar kind of transport barrier appears in the shearless region, one that is associated with a proper route of transition to chaos. These barriers have been identified in symplectic nontwist maps that model such zonal flows. We use the so-called standard nontwist map, a paradigmatic example of nontwist systems, to analyze the parameter dependence of the transport through a broken shearless barrier. On varying a proper control parameter, we identify the onset of structures with high stickiness that give rise to an effective barrier near the broken shearless curve. Moreover, we show how these stickiness structures, and the concomitant transport reduction in the shearless region, are determined by a homoclinic tangle of the remaining dominant twin island chains. We use the finite-time rotation number, a recently proposed diagnostic, to identify transport barriers that separate different regions of stickiness. The identified barriers are comparable to those obtained by using finite-time Lyapunov exponents.

Controlling chaos in wave-particle interactions.

We analyze the behavior of a relativistic particle moving under the influence of a uniform magnetic field and a stationary electrostatic wave. We work with a set of pulsed waves that allows us to obtain an exact map for the system. We also use a method of control for near-integrable Hamiltonians that consists of the addition of a small and simple control term to the system. This control term creates invariant tori in phase space that prevent chaos from spreading to large regions, making the controlled dynamics more regular. We show numerically that the control term just slightly modifies the system but is able to drastically reduce chaos with a low additional cost of energy. Moreover, we discuss how the control of chaos and the consequent recovery of regular trajectories in phase space are useful to improve regular particle acceleration.